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Wonderful life at weak Coulomb interaction: increasing of superconducting /superfluid transition temperatun 

Abstract. We have shown that in systems where the Coulomb interaction is strongly 
suppressed, the superconducting transition temperature may be enhanced by disorder 
close to the Anderson localization transition. This phenomenon is based on the 
enhancement by disorder of the matrix element of attraction in the Cooper channel. 
For enhancement to take place one needs (i) strong disorder which makes the single- 
particle wave functions strongly inhomogeneous in space and (ii) strong correlation 
of the patterns of inhomogeneity for different wavefunctions. One case where such 
correlation is well known is the system close to the Anderson transition. We review the 
notion of multifractality of wavefunctions in this region and show how the enhancement 
of Tc arises out of the multifractal correlations. 
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1. Introduction 

The statement which is known under the name of the "Anderson theorem" and belongs 
to P.W.Anderson [1] and Abrikosov and Gorkov [2], reads that the non-magnetic disorder 
does not change the superconducting transition temperature Tc. As a matter of fact 
this statement has a status of a theorem based on the normahzation and completeness 
of the set of single-particle wave functions only if the paring amplitude A does not have 
any variations in space or the latter can be averaged independently of the fluctuations of 
the single-particle wave functions. This is not true for sufficiently strong disorder where 
the joint effect of disorder and Coulomb interaction leads to suppression of Tc. This 
effect studied in detail by Finkelstein [3] was actually predicted earlier [4j as the leading 
correction to the mean field transition temperature in 2D disordered superconductors: 



where Agg = lyoV is the dimensionless constant of local interaction Vg = V, g is the 
dimensionless conductance, uq is the DoS at the Fermi level, D is the diffusion constant 
and uq is the Debye frequency. 

For the screened Coulomb interaction in the universal limit X^s = 1 and Eq.([T]) 
corresponds to suppression of Tc by disorder. 



However, if one assumes that the Coulomb interaction is absent (or strongly suppressed) 
and the only interaction that remains is the attraction in the Cooper channel 
characterizing by the small dimensionless constant Aeff = —A, the same Eq.([T]) together 
with the BCS relation A~^ = ln{uQ/Tc) results in: 



Thus in the absence of Coulomb interaction one obtains the enhancement (w.r.t. the 
BCS result) of Tc by disorder by the same token as the suppression in the universal limit 
of the screened Coulomb interaction known as the Finkelstein effect [3]. 

It was not understood in the eighties that the origin of both suppression and 
enhancement is the weak multifractality of single-particle wave functions in 2D metals. 
The main goal of this paper is to make this connection physically transparent. 

2. What is the multifractality? 

While the earlier period of Anderson localization was mostly devoted to study of the 
scaling behavior of the localization or correlation length as a function of proximity to 
the mobility edge E — Ec or to critical disorder W — Wc, the notion of multifractality 
of random critical single-particle wave function ?/'n(r) was introduced by Wegner [5] 




(1) 




(2) 




(3) 
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almost immediately after the formulation of the scaling theory of Anderson transition. 
The definition was given in terms of the inverse participation ratio and its moments: 

J d'r ^ L-^'-''"'\ (4) 

where (...) denotes disorder average, L is the size of the system and dq are the fractal 
dimensions. Eq. (jl]) covers all the three possible phases. The metal phase is characterized 
by the wave functions which occupy all the available volume with |-?/'„(r)p ~ L^'^ by 
normalization. In this case one immediately concludes from Eq. (j4]) that all dq are equal 
to the space dimensionality d. In the Anderson insulator |'?/'„(r)p is of the order of 
|'0„(r)p ~ ^^'^ with the probability {C,/LY for a point r to fall inside the localization 
volume ^"^j and is very small otherwise with the overwhelming probability 1 — {^/LY. 
Thus the leading term in Eq.dl]) at large L is of the order of L'^^~'''^{^/ LY ~ ^-(9-1)'^ 
and hence it is independent of L, which formally corresponds to all dq = 0. In the 
critical region close to the Anderson transition the set of critical exponents dq is the an 
important and non-trivial characteristic of this transition, in addition to the well known 
exponent u of the localization length ^ oc \ W — Wc\~'^. 

Another way to characterize multifractality is to define the spectrum of fractal 
dimensions /(a) so that the volume in space where the amplitude |'?/'n(r)P scales 
like L~" is Va (X L'f^"\ This approach explicitly takes into account the main feature 
of multifractality: the hierarchy of regions in space with different scaling of with 
the system size, each region being a fractal of the Hausdorf dimension /(a). With this 
definition we have for the integral of Eq. (jl]) : 

j^-[q-i)d, ^ da L''^''+f(''\ (5) 





By doing this integral in the saddle-point approximation valid at large In L one concludes 
that the quantity = {q — 1) dq and /(a) are related by the Legendre transformation: 

df{a) 

Tq = aq-f{a), —^ = q. (6) 

The validity of scaling Eq.Q at the Anderson transition point was demonstrated in 
numerous numerical studies and analytically by e = — 2 expansion [5] of the nonlinear 
sigma-model (weak multifractality (i — (i^ << c?) or by the virial expansion method [6] 
(strong multifractality dq « d). 

For further discussion it is very important to know the correlation of two single 
particle wavefunctions with the energies En = E and E^n = E + u which determines 
the matrix element of the local interaction. At the Anderson transition point this 
matrix element appears to be a power-law function of the energy difference u as was 
conjectured by Chalker |^, checked by direct diagonalization of the lattice Anderson 
models and critical random matrix ensembles [8] and finally proven analytically for the 
critical random matrix ensemble [9]: 

= d'r{\M^)\'\i^E+M)\') ~ [^y, (7) 
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where the exponent fi is closely related |7j with the fractal dimension d2'- 

The scaling Eq.(I7]),(|8]) in the frequency domain reflects the critical nature of single- 
particle states near the Anderson transition. However, it is valid well outside the 
mobility edge. In Fig.l it is shown the result of calculation of the correlation function 
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Figure 1. (color online) The correlation function Eq.([7]) as a function of cj from 
direct diagonalization of the 3D Anderson model with critical disorder W=16.5 for 
disorder strength W=2 (green), W=5 (red) and W=10 (blue), and different system 
sizes L=8(diamond),L=16 (circles) ,L=20 (squares). At = 10 and L — 16,20 the 
critical power-law is well seen for {i^o^'^)^^ — < uj < Eq. 



Eq.(I71) by direct diagonalization of the 3D Anderson model in the metallic phase with 
sub-critical disorder. One can see that even relatively far from the critical disorder 
(Wc = 16.5), at W = 10, the power-law critical behavior is well seen. It is saturated 
when the w-dependent "resolution length" = {uqu)~^^'^ becomes larger than either 
the correlation length ^ oc {Wc — W)~'^ or the system size L. Even at disorder as small as 
W = 5 the enhancement of given by the factor (Eq/u)^ in Eq.(I71) is quite significant. 
Only at W = 2 when ^ < i = (uqEq)^ it disappears, and we reach the limit L'^M^ = 1 
expected for the extended single particle wavefunctions which occupy all the available 
space. 

Such a behavior of the correlation function Eq. ([7]) suggests a stricture of the typical 
single-particle wavefunction at the sub-critical disorder which is sketched in Fig.2. It 
consists of the "fractal regions" of the size of the "correlation radius" ^ appearing in the 
scaling theory of localization. Inside such regions all the correlations are identical to 
those in the critical wavefunction. As the energy difference u decreases the "resolution 
length" = (uqu)'^ increases and at co < = (t'o^'^)"^ it exceeds the correlation 
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Figure 2. (color online) The cartoon of the amplitude landscape of the sub-critical 
single-particle wavefunction. The correlation length ^ oc {Wc — W)~'^ sets the typical 
size of the "fractal region" where the wavefunction correlations are the same as in the 
critical state. If the w-dependent "resolution length" = {i'QLd)~d is smaller than ^ 
the critical power law in Eq. ([T]) is well seen. For larger resolution length one observes a 
saturation in the w-dependence which is typical of metal. However, the fact that there 
are many peaks (dark regions on the plot)in the wavefunction landscape where the 
amplitude is pretty high, leads to the enhanced saturation level. The exception is the 
case of a very weak disorder when the correlation length ^ ^ £ is smaller than the size 
i = {i'oEo)~d of the "pixel" which is only slightly larger than the atomic scale. Only in 
this limit the Anderson theorem may apply. The similar cartoon holds for the localized 
states at super-critical disorder. In this case the amplitude is exponentially small 
outside the localization radius ^ but has the same fractal structure inside this radius. 
That allows to speak about the " multifractal metal" and " multifractal insulator" |8] . 

length ^. Only for such u the correlation function Eq.(I71) may show that the state is not 
critical but extended. Thus u = 6^ marks the onset of saturation in Fig.l for L 3> ^. 

It is important that at large u = Eq all curves in Fig.l merge, in agreement with 
Eq.([7]). This means that for u < Eq the matrix element is enhanced compared to 
the ideal case of totally extended wavefunctions which occupy all the available space 
without peaks and holes at a scale much larger than the Fermi wavelength. In order 
to better understand the origin of this enhancement we emphasize that it is caused by 
the inhomogeneity of the amplitude |'?/'(r)p in space. The presence of "holes" with very 
small amplitude in certain regions implies appearance of peaks in other regions, as the 
integral of the amplitude over entire system of the size L is equal to 1 by normalization. 
The key point is that the positions of the peaks in two different wave functions are highly 
correlated for the critical and slightly off-critical states. Contrary to naive expectation, 
even at large energy difference u which is only slightly smaller than the upper energy scale 
for fractality Eq (for the 3D Anderson model it is approximately 1/3 of the bandwidth). 
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the correlation is pretty high. This leads to the very slow decay of the matrix element 
~ uj~^ (yU < 1) with increasing u [TOl Hi]. Note that the ultimate reason for such a 
correlation is the presence of valleys in the random potential. All the wave functions are 
the solutions to the Schroedinger equation in the same random potential. So, the peaks 
know where they predominantly want to be situated. Yet, the criticality matters, as for 
localized states the similar correlation is absent, perhaps because there are too many 
possibilities to arrange a localized state of small localization radius. This is the reason 
why Eq.([7]) can be formally extended to ideal metal by setting d2 = d but it cannot be 
extended to hard insulator which formally corresponds to c/2 = 0. 

3. The minimal model in real and in Fock space 

As has been already mentioned in Introduction, in the situation when Coulomb 
interaction can be neglected the minimal model consists of the single-particle 
Hamiltonian with disorder and the local attractive interaction |12] : 

H=Y: er,r'^l(r)vI/,(rO-t/E4W^l(r)^i(r)^t(r). (9) 

r,r',a i" 

In the case where disorder is strong it is very useful to switch from the coordinate 
representation to the representation of the exact single-particle states where all effect 
of disorder is taken into account [131 [Hj. This is done as usual by substitution \1/Q,(r) = 
J2n V'nli") c„,a of the \l/-operator in terms of the exact single-particle wavefunctions V'nli") 
and the creation/annihilation operators c}^^ and Cn,a of electrons in the corresponding 
state n with the spin a. As the result the single-particle part of the Hamiltonian becomes 
diagonal but the interaction part produces a variety of terms 

H = ^En (^n,oPn,a - U ^ Mnkmlci^fcl^^Cm,lCLf, (10) 
",Q! n,k,m,l 

where 

M„,„z = J d^r^:(r)^^(r)^„(r)^Kr). (11) 

These terms may be divided into two major groups. The first one consists of the terms 
where all four indices are the same n = k = m = l or there are two pairs of equal indices. 
The examples are M„ = Mnnnn and M^m = Mnnmm, as well as Mnm = Mnmmn- They do 
not change the parity of occupation of a single particle state but can only move a singlet 
pair from state n to state m {Mnnmm) or to flip the spin of the singly occupied state 
{Mnmnm)- AH the othcr terms may convert the singly-occupied state into the empty or 
doubly occupied states and vise versa. 

In the region of localized single-particle states the largest matrix element is M„, 
all the other matrix elements are smaller because different single-particle states rarely 
overlap. When the matrix element — ?7M„ is large enough, it favors a local pair formation 
by producing a gap between the even many-body states in which all the single-particle 
states are occupied by a singlet pair or empty, and all other {odd) states. In this case 
one may neglect the presence of odd states whatsoever and introduce bosonic operators 
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^li ~ '^it'^lii- C)ne can easily check that = 0, so that we introduced a hard-core 

bosons. In addition to that, for even states = cj^^c^t ~ '^lii'^ni- 

The minimal model in the Fock space is obtained by neglecting of all terms except 
for those proportional to M„ and Mnm- The corresponding Hamiltonian takes the form: 

H = Y.2Enhihn-UY.Mnn,hihm + h.C. (12) 

n n,m 

It describes physics of local pairs with hard-core interaction attached to single-particle 
states n with random energies 2En and hopping due to the attractive interaction of 
original electrons. The term proportional to M„ is included into the chemical potential. 
Obviously, the hopping favors establishing a delocalized phase which at low enough 
temperature must be superfluid. Disorder, in contrast, is trying to localize the pairs. 

Eq. (fT2|l admits also a spin representation in which the hard-core nature of bosons 
is automatically included. To this end one introduces the spin-1/2 operators S:^ = bj^ = 
cl^cl^ and = b^bn - | = ^ [clt^nt + 4ii'^ni - !]• Then Eq.([T2]) takes the form: 

H = J22E^S:,-UY. M„„ S+S- + h.c. (13) 

n n,m 

In this language the Ising order corresponds to localized pairs and the XY order 
corresponds to a superfluid phase. 

So far we neglected the term proportional to Mnmmn which is also operating in 
the even sector of the Hilbert space. It can be easily taken into account by adding the 
following term to the Hamiltonian: 

-UY^Mnm S:,Sl, = -UY. Mnm bibJjOm + UV'^ ^ 6+ 6„. (14) 
n,m n,m n 

This term describes interaction (attraction) of bosons at different points in the Fock 
space. It favors phase separation with fixed pair occupation number 1 or and thus 
disfavors establishing a coherent superfluid phase where the occupation number strongly 
fluctuates. This has an effect of reducing the superconducting/superfluid transition 
temperature. Such a suppression of by a factor of order 1 [16] does not change, 
however, the main conclusions of this paper. 

Among other omitted terms the most significant are the following: 

-U J2 ^^nnrnl blCmiCi^ + h.C. (15) 
n,m,l 

which describe dissociation and creation of pairs out of singly occupied single-particle 
states. They can be safely neglected in the region of sufficiently strongly localized single- 
particle states but have to be taken into account in the region of critical and sub-critical 
states where they determine the Ginzburg number Gi. 

4. Mean field approximation in the Fock space 

The minimal model Eq.(IT3l) is a good starting point not only in the region of localized 
single-particle states {the "pseudo-gap" region) but also in the critical and sub-critical 
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region, although in this case its derivation from the minimal model in real space Eq.Q 
is not justified by a small parameter. Indeed, the standard mean-field treatment of the 
spin-model Eq.f llSp yields the following equation for the critical temperature [T5] : 

A r tanh (^] 
Ae = - dE' Ae Me-e', (16) 



2 J ^ E' 

where is given by Eq.Q and A = Uuq is the dimensionless attraction constant. 

At a weak disorder = 1 for |ci;| < oje, and zero otherwise, and one can 
immediately recognize the standard BCS equation for T^, albeit with T instead of 2T. 
This is the price of restriction to the odd sector of the Hilbert space. In case when 
there is no gap between the many-body states of the even and odd sectors, the latter is 
thermodynamically relevant even if dynamically both sectors are totally decoupled. We 
will show, however, that the accuracy of the mean-field approximation on the metallic 
side close to the Anderson localization transition is up to a constant pre-factor of order 
1. With this uncertainty the difference by a factor of 2 discussed above is beyond the 
accuracy of the mean-field approach. 

Let us apply the mean-field equation (fT6!) to compute the in the critical region. 
Plugging Eq.(I7j) in Eq.( fT6i) and making a rescaling E ^ TcE, E' T^e' one reduces 
the problem of to that of the eigenvalue A^^ = A^^ {Tc/Eq)^ ~ 1 of the dimensionless 
linear integral operator. Then we immediately conclude that: 

T.^EoXK /^ = 1-^- (17) 

d 

5. Ginzburg number and the accuracy of the mean-field approximation 

As is well known, the thermodynamic fiuctuations of phase and fiuctuation of the local 

Tc due to disorder restrict the region of validity of the mean-field approximation. To 

take into account correctly the first effect in the region of extended single-particle states 

one needs to include dissociation processes Eq.flT^. Without such processes, for the 

minimum model in the Fock space Eq.( !T3|) . the mean-field approximation is exact in 

this region, as the matrix element Mnm couples to an infinite number of states in the 

thermodynamic limit. When dissociation processes are properly accounted for [16] one 

arrives at a remarkable result: 
AT 

Gi = ^ ~ 1, (18) 

where ATc is the "fiuctuation region". The effects of multifractality cancel out in the 
Gi number, and it appears to be a universal number of order 1 close to the Anderson 
localization transition |14] . 

In the region of sufficiently strongly localized single-particle states, the Gi number 
is mainly controlled by the effective coordination number K of states coupled by the 
matrix element Mnm in the minimal spin model Eq. f|T3|) . As one goes away from the 
Anderson transition, this number increases. One can show [16] that in this region: 



= Gil + (19) 



Wonderful life at weak Coulomb interaction: increasing of superconducting / superfluid transition temperatun 



where Gii and Gi2 are universal numbers of order 1. This effect of decreasing of 
coordination number with increasing disorder invahdates the mean-field approximation 
at sufficiently strong disorder and finally leads to the superconductor to insulator 
transition. 

Note however, that in the vicinity of the Anderson transition Eq. (|T8l) holds true. 
Given the relationship between the Ginzburg number and the energy scale J associated 
with the phase rigidity: 

rn{MF) 

J = (20) 

Gi3 

one concludes that the scale J and the mean-field transition temperature Tj*^^) given by 
Eq. ffTTl) differ by a universal factor of order 1. This implies that the phase fluctuations 
may reduce the true superconducting transition temperature by at most a universal 
factor of order 1 compared to the mean- field result Eq.([T] 



6. Virial expansion method 

In order to describe the behavior of transition temperature in the region of localized 
single-particle states where the mean-field approach breaks down, we exploit [IE] the 
idea of virial expansion applied to the spin-model Eq. (fT3|) . To this end we analytically 
express the Cooper susceptibility Xc'{T) (the response of the spin to the infinitesimal 
perturbation A S^) as a sum of contributions of clusters of M coupled spins: 

oo 

Xc{T) = E Xm{T), (21) 

M+l 

The true transition temperature corresponds to temperature Tc when this series becomes 
divergent. According to D'Alembert criterion this happens when: 

hm ^M±iM = 1. (22) 

In view of rapid increase of complexity of the analytical expressions for xm with 
increasing M, we truncated the series Eq. fl21l) and used the operative definition of in 
terms of small-cluster contributions and Tjp'^: 

Xi(T«) = X2(T«), X2(Tf ) = X3(Tf ). (23) 

By evaluating xi, X2 and xs taking the energies En and matrix elements Mnm from exact 
diagonalization of the Anderson model on a 3D lattice and solving Eq.( |23l) numerically 
we were able to find statistics of TJ^^ and T^^^ at different positions of Fermi level relative 
to the mobility edge. The average TJ^^ and Tj^^ appear to be in good agreement with 
each other which gave us a reasonable confidence in the convergence of the procedure. 
In Fig.3, the results of virial expansion calculation for the average (TJ^^) are presented 
as a function of the Fermi level position. It is remarkable that the truncated virial 
expansion method is in good agreement with the results of mean-field calculations in the 
region of extended and weakly localized states and gives substantially smaller transition 
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Figure 3. (color onliiie)Transition temperature as a function of tlic Fermi energy 
position relative the mobility edge Ec for the fixed dimensionless constant A ~ 0.08 
(blue circles) and for the fixed attractive interaction U — 1.7i(green circles), where t 
is the hopping integral of the Anderson model. In the insert: energy dependence of 
the mean density of states. Large energy corresponds to region of localized single- 
particle states. Inside this region the mean- field (MF) transition temperature is 

(2) 

substantially larger than that (Tc = {Tc )) found from the truncated virial expansion 
series. The Cooper susceptibilities X2{T) and XsiT) were calculated numerically using 
exact eigenfunctions and eigenvalues of the Anderson model with strength of disorder 
W ~ 4, size of system L ~ 20"^ and the Debye frequency ujq = 0.5. 



temperature for stronger localized single-particle states. This proves that the virial 
expansion method, even with only 3 terms retained in the series Eq. fl22|) . captures the 
decrease of the effective coordination number K and is not equivalent to the mean-field 
approximation. 

7. Enhancement of Tc near the Anderson transition. 

Equation f|T7|) shows that when the Fermi energy is at the mobility edge E = the 
critical temperature behaves as a power law of the dimensionless attraction constant A. 
Numerical simulations on the 3D Anderson model of localization give (see Ref. [16] and 
references therein) the value of the fractal dimension d2 = 1.29 ±0.1. Thus the exponent 
in the power law dependence of T^, on A is close to l//i ~ 1.8. For small A one always has 
A^ *^ 3> exp(l/A). In addition to that the characteristic energy Eq is of electronic nature 
and it is typically higher than the Debye frequency uq. Therefore we conclude that 
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at small A the transition temperature given by Eq.f ll7p is much higher than the BCS 
transition temperature Tc ~ a;oexp(— 1/A). This enhancement of Tc is also obtained by 
a virial expansion method which is completely independent of the mean-field approach. 
It can be easily traced back to the enhancement of the matrix element Mnm shown in 
Fig.l. 

In order to figure out how general is the effect of enhancement of at negligible 
Coulomb interaction, let us consider the case of weak multifractality 1 — d2/d ^ 1. This 
case is relevant for the 2D metal [171 [18], where 

1-^ = -, 9 = 27x^1^,0. (24) 

2 g 

Expanding the matrix element Eq.^ in the mean-field equation Eq. f[TB]l up to the 
leading correction 'm. 1/g ^ 1 + ^ Iyi^Eq/uj) we obtain the 1/g correction to the 
transition temperature: 



5Z 



c 



2g \ TJ \uJoTj g KT^J ' ^'^^^ 

We obtained the correction of the same order as in Ref.||4j given by Eq.(j3]) thus proving 
that our analysis at A ^ 1 qualitatively applies also to the case of weak multifractality 
in 2D metals. Then Eq.([2]) tells us that the really crucial assumption for enhancement 
of Tc by disorder is the suppression of Coulomb interaction. Recently [19] this statement 
was confirmed by the RG analysis of the Finkelstein-like nonlinear sigma-model in d = 2 
and d = 2 + e where only the short-range interactions were taken into account. 

There are certain systems [20| |2T] where one may suspect such a suppression. 
Moreover, the enhancement of just before the onset of the insulating behavior was 
observed in one of them with the dependence of Tc on doping highly reminiscent of the 
Fig.3. More investigations are needed to confirm or reject this conjecture. In any case, 
the search for systems with suppressed Coulomb interaction (except for obvious case 
of cold neutral fermionic atoms with attractive interaction) is a challenging task which 
may open a door into a wonderful world without Coulomb interaction. 
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